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ABSTRACT 


A campaiison of data from all previous energy loss experiments per- 
formed at the Naval Postgraduate School with new theoretical predictions 
is presented in this thesis. In addition, the data have been extended 
to include beryllium. With this extension, experiments have now been 
conducted on materials ranging in atomic number from 4 to 82. Agreement. 
between experiment and theory is excellent for the most probable energy 
loss. However, theoretical values for the half-widths of the energy loss 
distributions generally are small compared to experiment for thicknesses 
less than 5 g/cm? and large for thicknesses greater than 5 g/cm*. These 
experiments were conducted in the energy range from 50 to 100 MeV, with 
Maes: thicknesses ean 0.7 to 5.9 g/em*. These effects were found to 


be independent of atomic number or incident energy. 
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I. INTRODUCTION 


Fast electrons which traverse a layer of material will exit with 
less energy than they possessed on entry. “he energy loss distribution 
can be characterized by two parameters, the most probable energy loss 
(Q,) and the half-width (HW), which is the full width at half maximum 
of the distrikution curve. These quantities are the basis for compari- 
son between experimental and theoretical results. Losses can be 
attributed to two major effects, ionization and radiation (bremsstrah- 
lung). Ionization losses for fast electrons have been calculated by 
Landau [1], Williams [2], and many others. Early Soecinanes resulted 
in good agreement with theory for the energy loss, but these results 
gave a half-width wnich was grcatox than theory at energies ahove 1 MeV, 
Blunck and Westphal [3] modified previous theoretical calculations by 
introducing radiation losses. This additional energy loss is super- 
imposed on the ionization loss and results in a larger half-width and 
enhancement of the tail of the distribution. Other losses have been 
discounted as unimportant relative to oneneen and radiation. All of 
the theoretical calculations are limited to “thin" absorbers, i.e., 
those in which 2, is small compared: to the incident electron energy. 

Several measurements of energy loss for high energy electrons through 
Various materials have been conducted at the Naval Postgraduate School. 
All were performed within an energy range of 50 to 100 MeV. These 
include the works of Bumiller, Buskirk, Dyer, and Miller on aluminum 


[4,5]; Goodwin on aluminum and copper [6]; DeLeuil, Raynis, et al on 








aluminum, copper, and lead [7,8]; and Mosbrooker and Sandquist on tin, 
lead, and gadolinium [9]. With the extension to beryllium, data are 
now available on various materials with atomic numbers (Z) between 4 
amc G2. 

Experimental consistency with past experiments was insured by 
rerunning several thicknesses of lead absorbers. Numerous and exten- 
Sive checks were made on the computer program (Appendix D) written by 
Mosbrooker and Sandquist [9] to insure that it was consistent with 
Blunck and Westphal [3] theory as modified by the beam folding tech- 
nique (Appendix C) which accounts for the energy spread of the incident 
beam. 

Comparison of all experimental data taken at the Naval Postgraduate 
School with the theoretical results of Blunck and Westphal [3] is pre- 
sented in this thesis. The data on beryliiium are entirely new; the 
older data have been presented before, but the theoretical values given 
here are new and represent calculations macle with the computer program 


of Appendix D. 








IT. THEORETICAL CONSIDERATIONS 


The Blunck and Westphal theory of the distribution for the energy 
loss of a beam of monoenergetic electrons in passing through a layer 
of absorbing material assumes that the energy loss Q is small compared 
to the initial beam energy, E.. Let W(Q)dQ be the probability of energy 
loss between Q and Q + dQ, and X be that portion of the loss Q due to 
radiation. Hence, the ionization loss is Q - X. Considering these two 


loss processes, the probability of energy loss is 


Q 
W(Q)dQ = 7 W, (Q-X)W, (x) dxdQ (1) 
O 


where Wy and We are the energy loss distributions for ionization and 
radiation respectively. 

The Landau equation [1], as modified by Blunck and Leisegang [10], 
is used for the energy loss distribution due to ionization. The distri- 


bution is expressed as a function of Landau's dimensionless parameter 


as follows: 


W.(0)dO = (1) dk = 2 oe a rl (2) 
ns ~ exp sb eee een as OE 0 
: n b? + y 2 b? + y.2 
a el 
_ 20 2 | 
where } = =~ + In | 1.116 (3) 


The terms used in equations (2) and (3) are defined as follows: 
Cir Ynr and A, are constants given by Blunck and Leisegang [10] and 
are used to fit Landau's distribution to a sum of gaussian functions. 


mews the absoroer thickness in an. 








The quantity "a" is a function of the atomic number Z, the atomic 
weight A, and the density p of the absorber; and 8(=v/c) of the 


electrons: 


MeV/cm. (4) 


The quantity b* is a correction to Landau's theory given by Blunck 


and Leisegang [10]: 





oe S20 . Tn io 2 a ; (5) 
aR Z 
m Lae) 


where the summation is over the ionization potentials of the atomic 
electrons, and NT is the number of electrons with ionization potential 
a2 e 
m 

Q is the average energy loss due to ionization (no radiation) for 


electrons of incident energy Eas and is given by Sternheimer [11,12,13] 


as follows: 


Te: m 
ae S 7 2 _ _ a S 
Q= “42 B+ 0.43 + Jn E. B c a. (xX, log, gp/mc) Mev (6) 
where t is the thickness in g/cm’, and the constants Aj, B, Cy Xr ag, 


and m. are parameters of the absorber material. These parameters for 
tin, lead, copper, beryllium, aluminum, and various other materials are 
listed in reference [13]. The parameters for gadolinium are not listed 
but were interpolated by Mosbrooker and Sanquist [9]. 

For We, the energy loss distribution due to radiation, Blunck and 


Westphal [3] give: 


a aR dQ 
We (Q)dQ = BaR (Q/E, ) 76) (7) 


10 








where 


Me eee 8 22 Pa \ eae ay cm7} (8) 
A 71/3 
ae: _ 1 
and B 1s a normalizing factor = ioR iD 


The distribution of total energy loss according to Blunck and 
Westphal is obtained by putting equations (7) and (2) into equation (1) 
and performing the required integration. The result is the energy loss 
Gistribution for a single electron of incident energy E.. For comparison 
of theoretical and experimental values, this distribution function was 
used, with corrections to account for the finite energy width of the 


incident electron beam. This treatment is described in Secticn IV and 


Appendix C. 
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TII. EXPERIMENTAL PROCEDURE 


Beams of electrons ranging in energy from 53 to 95 MeV were obtained 
by using the LINAC of the Naval Postgraduate School. Electrons were 
elastically scattered from the accelerator beam through 90° by a thin 
(1 mil) alumimm foil. The electrons passed through the varying 
absorber thicknesses and were energy analyzed by a 120° magnetic spec- 
trometer described by Kenaston, Luke and Sones [14]. 

The general experimental ,arrangement was similar to that used by 
Mosbrooker and Sandquist [9]. (See Figure 1). The energy of the inci- 
dent electron beam a determined by uSing a nuclear magnetic resonance 
probe to measure the field of the deflection magnet of the LINAC. The 
scattered electrons are countea on a ten channel piastic scintillator 
counting system that consists of ten aoe counters and a single backing 


counter which is operated in coincidence with the front counters. 
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Figure 1. Experimental Arrangement 
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Beryllium absorbers of thicknesses ranging from 0.738 g/cm? to 
5.908 g/an* were mounted on a remtely controlled ladder that could be 
moved up and down vertically in the electron beam. This allowed us to 
collect data cn one absorber thickness, then move the ladder and rm a 
different thickness without having to enter the LINAC end station or 
change the LINAC beam in any way. The absorbers were placed at about 
Ewen from the scattering foil. 

A lead absorber of thickness 2.825 g/cm* was used to compare our 
experimental data with the results of Mosbrooker and Sandquist [9]. 
The results were in agreement within experimental accuracy. 

Since the electron beam scattered into the absorbing material was 
not monoenergetic, the energy distribution was measured before and 
after passing through the absorbers. The zero thickness absorber dis- 
tributions found by passing the beam only through the scattering foil 
had half widths ranging from 0.21 MeV to 0.46 MeV, or an energy reso- 
lution of about 0.45%. Following the measurement of these "zero peaks", 
the varying alosorber thicknesses were moved into place and an energy 
distribution was determined for each absorber thickness. 

Because of the limited quantity of beryllium available and space 
requirements in the target chamber, it was not possible to place all 
eight absorber thicknesses on one ladder. Therefore, each time the 
beam energy was changed or the ladder arrangement was changed, a new 
zero peak was measured. This accounted for the possibility that the 
Character of the beam might change. 

The data represent the number of scattered electrons detected by 
the coincidence counting system at the exit of the spectrometer as a 


function of the energy selected by the spectrometer. 


13 








Normalization was controlled by integrating the current fran a 
Secondary Emission Monitor. This instrument was mounted both up and 
downstream of the targets. There proved to be no significant difference 
in the width of the zero peaks due to the position of the monitor. Each 
data set (ten points) represents a given integrated charge, i.e., a 


certain number of electrons passing through the scattering foil. 
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IV. TREATMENT OF DATA 


The data for the experiment were taken with the ten channel counting 
system as described in Section III. To reduce the data, account must be 
taken of the characteristics of the counting system. These include the 
energy spread of the counters, which is about 0.3% between counters or 
3% for the entire system. Accordingly, the energy seen by the individual 
counters must be determined and also correction must be made for the 
differing efficiencies of the scintillators and associated electronics. 

A computer program has been written which determines the correct relative 
counter energies ate also corrects for counter efficiencies. Background 
corrections, which were quite small, were determined by periodically 
removing the scatterer ana apsorbers fran the beam. 

The reduced data were then plete. end from these plots the experi- 
mental most probable energy losses and distribution half-widths were 
measured. These results and their uncertainties can be found in Table I 
in Appendix A. | 

Since the scattered electron beam is not monoenergetic and the Blunc* 
and Westphal theory is. based on a monoenergetic incident beam, to compare 
theory and experiment the theoretical model mst be modified. This modi- 
fication is accomplished by representing the incident distribution as a 
histogram, each bin of which has a definite energy, and then unfolding 
this distribution into the theoretical calculations carried out by an 
IBM 360/67 computer. The unfolding technique is described in Appendix C 


and the computer program used is in Apperdix D. 
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The camputer program described in Appendix D also calculates the 
theoretical half-widths and most probable erergy losses. The program 
nomalizes the theoretical curve to the experimental data and plots 
both the theoretical curve and experimental data. These plots for 
beryllium are shown in Figures [3] through [26] in Appendix B. The 
measurable parameters were compared and the final results are tabulated 
in Appendix A. 

The camparison of all data from all previous energy loss experi- 
ments performed at the Naval Postgraduate School LINAC was accomplished 
by normalizing each set of experimental data to the Blunck and Westphal 
theory as calculated by the ccmputer program described in Appendix D. 
The comparison of the theory parameters so obtained to the experimental 


data can be found in Tables II through VI in Appendix A. 
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V. - RESULTS AND CONCLUSIONS 


The theoretical predictions and experimental results for the most 
probable energy losses and half-widths for beryllium can be found in 
Table I, Appendix A. Our theoretical predictions and the experimental 
results from previous theses completed at the Naval Postgraduate School 
for aluminum, copper, tin, gadolinium and lead can be found in 
Tables II through VI respectively in Appendix A. The theoretical data 
shown are for the Blunck and Westphal theory as modified by the beam 
folding technique. In all tables BW refers to Blunck and Westphal theory, 
2, represents the most probable energy loss, and HW is the distribution 
half-width. The colum headed % represents the percentage by which theoxy 
Be ceads experiment. 

Agreement between the Blunck and Westphal theory and our experi- 
mental results is excellent for the most probable energy loss. The 
average difference between experiment and theory is less than 0.5%. 
However, agreement for half-width is good only for thicknesses less than 
1 g/an? and is fair for thicknesses between 4 to 6 g/cm’. This result 
in half-width comparison is in contrast to previous results where it 
was concluded that the Blunck and Westphal theory predicts satisfactorily 
for thin targets where "thin" was defined as t (2) 173<13.5, where t is 
the absorber thickness (g/cam¢) and z is the absorber atanic number. 

The observed trend for HW (experimental) divided by HW (theoretical) 


as a function of t (absorber thickness) is plotted in Figure 2. 
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As can be observed from these plots the theoretical half-width 
tends to be smaller than experiment until one gets to thick absorbers 
(approximately 5 g/cm’) at which point it gets larger than experiment. 

We have considered two possible effects which might explain why 
this crossing of experiment and theory takes place. One effect is that 
this trend is caused by an assumption that is made in the theory. This 
assumption is that of constant incident enexgy (E,) , 1.e., that the 
energy loss is small campared to the incident energy of the electrons. 
However, in thick absorbers, this assumption is not valid, therefore 
the theory tends to overestimate the energy loss due to radiation 
(bremsstrahlung). This overestimate will cause an increase in the dis- 
tribution tail and hence increase the distribution half-width. ‘This 
trend in thick absorbers is not surprising in that the theory was 
OLiginally designea for thin absorbers. The fact that the agreement is 
excellent for 2, while not for half-width is also not surprising since 
an overestimate of this type would not affect the most probable energy 
loss significantly. Although this theoretical inaccuracy is undoubtedly 
a cause of some discrepancy in agreement with experiment, it is doubtful 
that it is the cause of all the disagreement. | 

The second possibility which we have considered affects the experi- 
mental results. Due to multiple scattering in thick absorbers the possi- 
bility existed that the size of the spectrometer opening was excluding 
the low energy portion of the distribution. However, calculations of 
the projected mean square scattering angle as outlined in Jackson [15] 
showed this angle to be approximately 0.30° while the spectrometer 
opening was approximately 0.90°. Therefore, it wae determined that 
Multiple scattering was not a factor in the disagreement of theory and 


experiment. 
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Discrepancies in the theoretical predicticns fran our program and 
previous results of Goodwin [6] and DeLeuil and Raynis [8] were noted. 
An attempt was made to resolve these differences. The parameters from 
Blunck and Westphal [3] were inserted into our computer program and 
without the unfolding technique described in Appendix C we were able to 
duplicate closely the Blunck and Westphal cistributions. On this basis, 
it is assumed that our theoretical results are the more accurate. 

Contrary to previous theses, the thickness at which theory and 
experiment began to exhibit a noticeable discrepancy was found not to 
be dependent on the atcmic number of the material. 

It should be mentioned that no attempt was made in this experiment 
to observe the very low energy tail of beryllium (approximately at 15 
MeV) where in previous experiments secondary effects had caused a peak 
memine assorber distribaticn. 

In light of what has been found in this thesis, further work is 
justified in several areas. First of all, the low energy tails of the 
various absorber distributions should be investigated to find out if 
secondary effects do cause a low energy peak in some distributions. 
Secondly, an effort should be made to improve the efficiency of the 
current computer program. Finally, an attempt should be made to extend 
the theory to thick absorbers where the energy loss is not small with 
respect to the incident energy, and also to include the effects of 


secondary processes. 
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APPENDIX A — TABLES 


Energy Loss Distribution Characteristics of Berylliun 
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Energy Loss Distribution Characteristics of Copper 
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Table IV. Energy Loss Distribution Characteristics of Tin 
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Table V. Energy Loss Distribution Characteristics of Gadolinium 
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Table VI. Energy Loss Distribution Characteristics of Lead 
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APPENDIX C —- BEAM FOLDING TECHNIQUE 


The theory of Blunck and Westphal [3] assumes that the beam of 
electrons striking the absorber is monoenergetic. The beam of elec- 
trons produced by the Naval Postgraduate School LINAC, or any other linear 
accelerator, for that matter, has a finite energy spread about a maximum 
or most probable energy point. The fact that monoenergetic electrons 
are not available to strike the absorber must be taken into account in 
computing theoretical predictions if a meaningful comparisam™m with 
experimental results is to be made. This has been done here by a 
technique termed "beam folding". | 

Beam folding is accomplished by a number of well defined steps. 
Energy distribution curves may be anproximated by histograms. These 
histograms consist of a series of "bins" of area W(Q) AQ, where Q is 
the energy loss and W(Q)AQ is the probability of loss between Q and 
Q+ AQ. Thus, each bin has an "address", Q, on an energy coordinate 
scale. To accomplish beam folding, the beam distribution must be known. 
This is observed experimentally and aeceecive teed in the computer by a 
gaussian curve of appropriate half-width. The energy at which the 
maximum occurs is established by the energy of the beam and the magni- 
tude of energy loss incurred by elastic scattering of the beam as it 
impinges on the thin aluminum scattering foil prior to striking the 
absorber. The width AQ is then selected for the predicted distribution. 
This must be a small, but finite number where numerical techniques are 
to be used. This same width is used to break the beam distribution into 
a histogram. This is illustrated in Figure 30(a). The reason the sare 


width is used for both distributions is to facilitate computer programming. 
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Each bin of the beam distribution is now treated as a monoenergetic 
beam with energy commensurate with its center location on the energy 
scale. Each bin also has a definite sears or weight, with the magni- 
tude of the center bin being unity. The formulae of Blunck and Westphal 
is now applied to each of these beams and an absorber distribution curve 
results for each, with a maximum amplitude proportional to the height of 
the appropriate beam distribution Waves Bach of the absorber dis- 
tribution curves thus obtained may be thought of as being plotted and 
added to previously determined curves, using the energy of the electrons 
aS a correlation point, as Detoues in Figures 30(b), (c) and (dad). This 
is accomplished on the eure by adding the contents of each bin of the 
same address and plotting the cumulative total. The bin of address E,. 1s 
deoicted in Figures 30(b), (c) and (d). Note that the beam distribution 
for histogram E4 is centered over the maximum point of the distribution. 
This is important as a false picture could easily be presented if the 
histogram were not symmetrical as the distribution would then appear 


skewed. 
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APPENDIX D - COMPUTER PROGRAM 
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